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^ ■ We investigate the average velocity of Brownian particles driven by a constant exter- 

^! nal force when constrained to move in two-dimensional, weakly-corrugated channels. 

^ I We consider both the geometric confinement of the particles between solid walls as 

'^ ■ well as the soft confinement induced by a periodic potential. Using perturbation 

i__il methods we show that the leading order correction to the marginal probability dis- 

Q ' tribution of particles in the case of soft confinement is equal to that obtained in the 

C/3 ! 

^ , case of geometric confinement, provided that the (configuration) integral over the 

c^ ■ 

^ I cross-section of the confining potential is equal to the width of the solid channel. We 

'^ ■ then calculate the probability distribution and average velocity in the case of a sinu- 

O ! soidal variation in the width of the channels. The reduction on the average velocity 

is larger in the case of soft channels at small Peclet numbers and for relatively narrow 

^ . channels and the opposite is true at large Peclet numbers and for wide channels. In 

Q>^ I the limit of large Peclet numbers the convergence to bulk velocity is faster in the case 

\0 ■ 

^ ' of soft channels. The leading order correction to the average velocity and marginal 

^ ! probability distribution agree well with Brownian Dynamics simulations for the two 

^ I types of confinement and over a wide range of Peclet numbers. 

>'■ 

X- 



I. INTRODUCTION 

The transport of Brownian particles confined to a channel with periodically varying cross 
section is relevant to a broad range of problems, from tracer dispersion in porous media 
to diffusion across entropy barriers in biological systems.-'^ In addition, the transport of 
suspended species under geometric confinement is ubiquitous in the rapidly growing field 
of microfluidics. Various separation microdevices, for example, are based on the effect that 
heterogeneous microstructures and geometric confinement have on the average velocity of 
suspended particles.- 

A number of studies have investigated the case of narrow channels, in which the charac- 
teristic length scale of the cross-section is much smaller than the length of a single period 
along the channel. In the schematic of a two-dimensional wavy-wall channel shown in Fig. 
[1] a narrow channel would correspond to e ^ 1. In a narrow channel, the description of 
diffusive transport in the absence of an external force can be simplified by reducing the 
dimensionality of the problem via the Fick- Jacobs approximation, in which the motion in 
the cross section is transformed into an entropic barrier to longitudinal transport.-"^"- In 
the presence of an external field a relatively simple extension of the Fick- Jacobs approxima- 
tion has been used to calculate the average velocity of Brownian particles.-i^i^ Alternatively, 
Laachi et al}^ used the standard long- wave asymptotic perturbation analysis and obtained 
analogous results for the average velocity to leading order in e. In fact, we have shown that 
the leading order term in a perturbation analysis is equivalent to a direct extension of the 
Fick- Jacobs approximation to the case of biased transport. ^^ In addition, and motivated by 
recently proposed partition-induced separation devices in which suspended species are trans- 
ported above a patterned surface,— we also considered the transport of Brownian particles 
confined by a relatively narrow channel in the potential energy landscape.— Remarkably, 
we found that the confining potential has identical effects on the transport of Brownian 
particles to those induced by a solid channel, to leading order in the aspect ratio e. 

In this work, we study the transport of Brownian particles through a two-dimensional 
wavy-wall channel in the limit of small variations in the width of the channel, that is A ^ 1 
in the schematic shown in Fig. [H On the other hand, we shall assume that, in general, 
e ~ 0(1). Following our previous work we consider both the geometric confinement by solid 
walls as well as the soft confinement induced by a periodic energy landscape. This means 
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FIG. 1. Schematic diagram of the channel. The upper and lower "boundaries" of the channel corre- 
spond to the solid walls in the case of geometric confinement or to two representative equipotential 
lines in the case of soft confinement. 

that the channel boundaries shown in Fig. [1] should be interpreted either as the solid walls 
or as two representative equipotential lines. In both cases we investigate the effect that a 
small perturbation on an otherwise flat channel has on the average velocity of the particles. 

II. MODEL DESCRIPTION AND PREVIOUS RESULTS 



Consider the transport of Brownian particles in the two different types of wavy channels 
discussed above. In one case, the particles are confined between two solid walls. In the 
other case, the particles move in a potential landscape that confines them in the direction 
perpendicular to the plane of motion. Specifically, in the case of geometric confinement 
we consider a two-dimensional, symmetric channel with half-width h{x, z) = h{x) in the 
z-direction and periodic in the x-direction, h{x + L) = h{x). In the case of soft confinement 
we consider a periodic, two-dimensional potential, V{x + L,z) = V{x, z), that confines the 
particles in the z-direction, i. e. V{x, z) — )■ +00 for z — ?■ ±00. In both cases, in the limit of 
negligible inertia, the motion of Brownian particles can be described by the Smoluchowski 
equation for the probability density P{x,z,t), 



dP 

'dt 



+ V ■:i = 5{x,z)5{t). 



(1) 



The probability flux, J{x, z,t), is given by 

_ 1 / — dV— dP\ - 1 / dV— dP\ - 

rj \ ox ox J 7] \ oz oz ) 

where F is a uniform external force driving the particles in the x-direction, 77 is the viscous 

friction coefficient, k^ is the Boltzmann constant, T is the absolute temperature, and we 

have used the Stokes-Einstein equation to write the diffusion coefficient as a function of 77, 

D = ksT/rj. Let us mention that in the case of geometric confinement we shall assume that 

there is no potential field present in the system, that is V{x, z) = 0. 

In order to obtain the asymptotic distribution of particles, relative to a single period of 

the channel, we first introduce the reduced probability density and the reduced probability 

for a more detailed discussion and Ref. 



current (see Refs. 



13 and 



14 



15 



for an analogous 



approach based in macrotransport theory), 

+00 
P{x,z,t)= Y, P{x + n,L,z,t), (3) 

+00 

J{x,z,t)= 2_] Jix + nxL,z,t). (4) 



The reduced probability is then obtained by solving Eq. ([T]) with periodic boundary 
conditions in x, P{x, z, t) = P{x + L, z, t), the no-fiux condition in the z-direction, and the 
normalization condition over a unit cell. The no-fiux condition imposed by the confinement 
depends on the type of channel. In the case of geometric confinement, the condition at the 
solid walls is 

3-n = r h'{x) ±r = (5) 

where n = h'{x)i ± fc is a vector normal to the channel top and bottom walls defined by 
z = =p/i(x), respectively. In the case of soft confinement, the no-fiux condition at the solid 
boundary is replaced by a far-field condition in J^, i.e. a vanishingly small probability 
density (and fiux) due to the confining potential, 

1 / dV— dP\ 

7] \ oz OZ J z^±oo 

In what follows we simplify the analysis by using dimensionless variables, with L as the 
characteristic length along the channel and, to be consistent with our previous analysis in 
the case of narrow channels, we choose eL as the characteristic length scale in the cross- 
section, where in this case e is not necessarily small. Then, we introduce the dimensionless 



variables x = x/L and z = z/{eL), as well as the dimensionless potential, V = V/{kBT), 
and probability density P{x, z, t) = eL'^P{x, z, t). The governing equation for the asymptotic 
probability distribution, Poo{x,z), becomes 
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0, 



(7) 



where the Peclet number is defined as Pe = FL/ksT. The periodic boundary condition in 



dimensionless form is Pooix + l,z) 



Poo{x,z), and the normalization condition becomes 



J rv^ f-t/iX' H/ lO -L . 



(8) 



where the domain Q is given by {(x,^;) : < x < 1, — oo < z < oo} for the case of 
soft confinement, or {{x,z) : < x < l,—h{x) < z < h{x)} for the case of geometric 
confinement by solid walls. Finally, the dimensionless no-fiux condition in the case of soft 
confinement is essentially the same as in Eq. ([6]), that is J^ — )■ for z — > ±oo. On the 
other hand, the no-fiux condition at the solid boundary becomes 



J . n = e'r h'(x) ± J 



0, 



(9) 



due to the different characteristic scales used to nondimensionalize x and z. 

After the asymptotic solution for the reduced probability distribution is determined we 
can obtain the average velocity along the channel by applying macrotransport theory,— 
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JZ-, dx dz 



Pe-^Voo-^ 
dx J dx 



dxdz. 



(10) 



Alternatively, the motion of Brownian particles in a viscous solvent in the limit of van- 
ishingly small inertia can also be described by the overdamped Langevin equations, 



dx 



dV 
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and 



dz 



dV 



V^ = -^ + V^iknTUt), 



(11) 



(12) 



where C{t) is a zero-mean, Gaussian white noise, with independent components satisfying 
the fiuctuation-dissipation theorem, (Cj(^i)Cj(^2)) = 2(5jj5(ti — t2)- In this case, the aver- 
age velocity along the channel can be evaluated from an ensemble average of independent 
trajectories. 



(v) 



lim 

t—^oo 



(x(t)) 



(13) 



A. Asymptotic analysis in the narrow channel approximation 

In Ref. lUJ 



we showed that, in the case of narrow channels (e ^ 1), the leading order 
solution of the reduced probability is equivalent to the probability distribution that is ob- 
tained from the extension of the Fick- Jacobs approximation to the case of biased transport. 
Moreover, we showed that the geometric confinement between solid walls has the same effect 
on the transport properties of Brownian particles as the confinement by a soft potential, to 
leading order in e. Specifically, the governing equations for both geometric and soft confine- 
ment cases are the same, provided that the channel width, w{x) = 2h{x), is equal to the 
following integral over the cross section, 

/■oo 

I{x) = / e-^("'")dz. (14) 



oo 



Note that, in the original work by Zwanzig,^ this integral was associated with a free energy 
A{x) given by I{x) = exp(— ^(2;)). It can also be considered as a configuration integral over 
all possible states (vertical positions of the particle) at a given position along the channel. 
The case of solid walls can be modeled as an infinite square-well, in which the potential is 
zero inside the channel and infinite outside. In this case it is immediate that I{x) = w{x) 
and the free energy is entropic. 

If the condition described above is satisfied the leading order correction to the average 
velocity is also the same in both cases. 



{vfj) = (1 - e 



-Pc\ 



(15) 



III. ASYMPTOTIC ANALYSIS IN THE SMALL PERTURBATION 
APPROXIMATION 

A. Soft confinement by a periodic potential 

We consider the transport of Brownian particles confined to a channel by a potential of 
the form 

V{x, z) = Voiz){l + XV^{x) + AV2(x) + ■■■)• (16) 

The potential Vo(^), by itself, would confine the particles to a straight channel, assuming that 
Vo{z) — )► 00 for z — )► ±00. The periodic functions Vn{x) (for n = 1, 2, 3 ■ ■ ■ ) therefore act as a 



perturbation to the channel width, where the magnitude of the perturbation is determined 
by A. Here, we analyze the case of small perturbations, that is A ^ 1. Therefore, we look 
for a solution to Eq. ([7]) in the form of a regular perturbation expansion for the asymptotic 
probability distribution, 

Poo{x, z) ~ Vq{x, z) + Api(a;, z) + AV(a;, z) A , (17) 

where each function "PrSs^^ V) (n = 0, 1, 2, 3 ■ ■ • ) is periodic in x, p^ [x = 0, z) = p^ [x = l,z). 
Analogously, we write a regular perturbation expansion for the probability flux, 

Joo(a;, z) ~ Jo(a;, z) + AJi(x, z) + A^J2(a;, z) + ■ ■ ■ . (18) 



Let us also write the perturbation expansion of the configuration integral in Eq. (iMj) ex- 
plicitly as 

I{x) = Jo (1 + A/i + 0(A2)) 

_ (19) 

e-^^'^'Uz (1 - Ayi(x)ro + 0{X^)) , 

where 

J —oo 

Substituting the perturbation expansion for the probability density into the governing 
equation, Eq. (^^, and equating like powers of A we obtain a hierarchy of coupled equations 
for the different functions Pn(x, z). The equations governing the basic solution, po{x, z), and 
the perturbation functions pi{x, z) and P2{x, z) are: 

e^— (ve - — ^ + — (-'^ - ^^ = (21) 

dx \ dx ) dz \ dz dz ) 

and 
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dx 



fdV2 dV, \ dp, 

Pe P2-Vg\ -1— Po + -1— Pi - 7^ 
\dx dx J ox 



(23) 



which correspond to the 0(1), 0(A), and O(A^) terms, respectively. Analogously, replacing 
the regular expansion for the flux in Eq. (jH]) we obtain, 

J^(a;,z ^±oo) = 0, for n = 0,1,2- ■■. (24) 



The normalization condition for the probabihty density corresponds to a basic solution 
Po{x, z) that is normalized to unity, 

(Po) = dx dz po = l, (25) 

Jo J-oo 

and to the perturbation functions having zero average (pn) = for n = 1, 2, 3 ■ ■ ■ . 

In order to calculate the average velocity, we first simplify Eq. fllOl) using the normaliza- 
tion condition and the periodicity of the probability distribution, 

(v) = Pe- dx dz J- Poo- (26) 

Jo J-oo C'X 

Replacing the expansion for the probability distribution into the equation above we obtain 
the general expansion for the average velocity to order A^, 

vo = Pe, (27) 

vi = - dx dz Vq—^ Po, (28) 

Jo J-oo "^ 

V2 = - r dx r dz Vo (^po + ^Pi] . (29) 

Jo J-oo \dx dx J 



The basic solution is completely determined by Eq. fl2T|) and the normalization condition, 
and we obtain: 

Po = |e-^«W, (30) 

-'0 

where Iq was defined in Eq. ( IT9l) . Clearly, this basic solution is independent of x and corre- 
sponds to the confinement by a straight channel in the absence of perturbations. Therefore, 
it is easy to see from Eq. ( 128|) that f i = due to the periodicity of the perturbation po- 
tential V^i(x). Then, in order to obtain the leading order correction to the average velocity, 
we first need to determine the leading order correction to the probability density, that is 
Pi{x,z). The eigenvalue problem associated with Eq. fl22l) and the representation of the 
general solution in terms of eigenfunctions are discussed in Appendix El and, in the next 
section, we will present the exact solution for the case of a parabolic potential confining the 
particles in the 2;-direction. Here, we limit the general discussion to the simpler problem of 
finding the correction to the marginal probability density at 0(A), that is 

/■oo 

Pi{x) = / pi{x,z)dz. (31) 



The governing equation for pi{x, z) is obtained by integrating Eq. (122|) over the cross section, 
which yields 

ax \ ax ax J 

where Ii{x) = —Vi{x)Vo was defined in Eq. ( 1T^ . Integrating Eq. ([32]), taking into account 
that f 1 = 0, and enforcing the normahzation condition we obtain 



Pii^) = -:^^:—r / e"^'^ -J7 ^C (33) 



e-P-^ - 1 A ^C 

i. Particles confined by a parabolic potential 

In order to explicitly calculate the leading order correction to the average velocity, we 

consider the transport of Brownian particles confined to a channel by a parabolic potential 

of periodically fiuctuating width. Specifically, the potential in non-dimensional variables is 

given by 

1 / z 



where 



(5 (x) = W ^ Q + A sin(27rx) j , (35) 

The equilibrium probability is given by a Gaussian distribution with variance (j{x) = 6{x). 
Let us mention that the factor \/2/y/TT included in the definition of 5{x) yields I{x) = 
1 + 2A sin(27rx), which will be convenient when comparing the results with the case of solid 
confinement. Clearly, in this case Iq = 1 and /i(x) = 2sin(27rx). 
The expansion of this potential is, 

V{x, z) = nz'^ (l - 4A sin(27rx) + 12A2 sin2(27rx) + ••■), (36) 

and the basic solution is the Boltzmann distribution with Vo{z) = vrz^ as the potential 
energy, 

p,{z)=e~''^\ (37) 



Therefore, the corresponding leading order marginal probability distribution is uniform. The 
first correction to the marginal probability is obtained by direct integration of Eq. fl33l) . 



Pi{x) = Att 



27r sin 2TTX — Pe cos 2tix 
47r2 + Pe^ 



(38) 



We can also write this marginal distribution in terms of a phase-shift 0o with respect to the 
confining potential, 

/Pe\ 
Pi{x) = 2cos0osin(27rx — (J)q); 0o = arctan I 7^ 1 • (39) 

Fig. [2] shows the marginal probability for different Peclet numbers. Interestingly, in the 
limit of small Peclet numbers we obtain pi{x) — ?■ 2 sin27ra;. This is consistent with a uniform 
distribution across a solid channel with a width given by w{x) = I{x) = 1 + 2Asin27rx. On 
the other hand, at large Peclet numbers, the distribution of particles integrated over the 
cross-section becomes completely out of phase with respect to the variations of the channel 
width, that is 0o — ^ 7r/2. However, in this limit the magnitude of the correction to the basic 
solution vanishes. 

As we discussed in the previous section, in order to obtain the leading order correction to 
the average velocity, V2, we need to calculate pi{x, z). The governing equation for pi{x, z) is 
obtained from Eq. (!22l) by substituting the expressions for the confining potential, Vq = vr^;^, 
and the leading order perturbation, Vi = — 4sin(27rx), 

e'^ (pep, - ^"j + ^ f-27Tzp, - ^) = [{lQe\' + 167r')z' - Svr] posin(2vrx). (40) 



D. -0.5 




FIG. 2. Marginal probability density pi{x) in the case of soft confinement by a parabolic potential, 
for different Peclet numbers. 

In Appendix |X] we show that, in the case of a parabolic potential, pi{x, z) can be written 
in terms oi po{z)Ho{y/nz) and po{z)H2{y/nz), where Hm{x) are the Hermite polynomials, 



10 



Hq(x) = 1 and H2{x) = (4x^ — 2). Let us note that the latter term, 2po{z)(27iz'^ — 1), 
integrates to zero over the cross section and thus it does not contribute to the marginal 
probability. Therefore, we propose the following solution, 

Pi{x,z) = [2cos(0o) sin(27rx — 0o) + i2nz'^ — l)(asin(27rx) + 6cos(27rx))] Poiz), (41) 

where a and b are undetermined constants. The proposed solution is clearly periodic and 
the normalization condition is also satisfied for arbitrary values of a and b, by construction. 
These two constants are determined by substituting Eq. (14T!) into Eq. (140|) and equating 
the coefficients of each of the orthogonal functions to zero. The final solution is, 

Pi(x, z) = 2 [cos(0o) sin(27ra; — 0o) + cos(0i)(27r2;^ — 1) sin(27rx — 0i)] poiz), (42) 



where the phase shift (pi is given by, 

01 = arctan 



e^Pe 



2(7re2 + r 



(43) 



Finally, the leading order correction to the average velocity, V2, is obtained from Eq. ( l29ll 

dz (vTj;^) (—8n cos(27rx)) pi{x, z) 



V2 



1 fOO 

dx 

-J —CO 

— 47r [cos(0o) sin(0o) + 2 cos(0i) sin(0i)] 

. p vr 2eHne' + 1] 

— avrre 



(44) 



_47r2 + Pe2 e4Pe2 + 4(7re2 + l)2_ 

Alternatively, we can calculate the effective mobility, defined as the ratio of the average 
velocity to the applied force. In dimensionless form, the effective mobility normalized by 
the mobility in bulk, /io = 1]'^, is given by 

2e\7Te^ + 1) 



Ho Pe 



TT 



+ 



+ 0{X' 



(45) 



_47r2 + Pe2 e4Pe2 + 4(7re2 + l)2_ 

Finally, and for comparison with the case of geometric confinement, we calculate the limiting 
behavior at small Peclet numbers. 



limM = i_2A^ 

Pe-s.0 Pe 



1 + 



27re2 



7re2 + l 



(46) 



and the asymptotic behavior at large Peclet numbers, 

Pe e^Pe^ ^ ^ 



(47) 
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B. Geometric confinement by solid walls 

Let us consider now the case of a solid channel that is periodic in x and symmetric about 
the xy-p\a.ne. The upper wall of the channel is given hj z = h{x), with 

h{x) = ^ + Xg{x), (48) 

where A is again assumed to be a small parameter. Note that this choice of channel width 
corresponds to Jq = 1 in the case of soft-confinement. The main objective is to determine 
the effect that such small oscillations in the width of channel have on the average velocity 
of the particle. Therefore, we seek a solution to the asymptotic distribution in the form of 
a regular asymptotic expansion in A, 

Poo{x, y) ~ Pq{x, z) + Api(x, z) + \^p2{x, z)^ . (49) 

Replacing the expansion above into the governing equation for the asymptotic probability, 
Eq. (jTj), and equating like powers of A, it becomes clear that in the absence of an external 
potential the governing equations for successive powers of A are all the identical. Specifically, 
the general equation for p„ is 

e^— fpe - — ^ - ^^ = (50) 

dx \ " dx J dz"^ 

Although the previous equation seems to indicate that the different functions pn{x,z) are 
independent, they are actually coupled to each other through the boundary and normaliza- 
tion conditions. Therefore, we first use the domain perturbation method^i^^i to transform 
the no- flux condition at the wavy walls, given by Eq. (|9]), to an asymptotically equivalent 
boundary condition at the two flat planes given hj z = ±1/2. Analogously, Eq. (|8]) is 
transformed into an equivalent normalization condition using a Taylor series expansion for 
the domain Q (a more detailed discussion is presented in Appendix |B]) . Finally, applying 
the domain perturbation method to Eq. (ITUl) we shall obtain the regular expansion for the 
average velocity of the particles. 

Let us first rewrite Eq. (fTO!) replacing Poo{x,z) by its regular expansion. 

Before we expand the domain of integration in powers of A, we note that the basic solution 
for the asymptotic distribution of particles is uniform, that is po = 1- This is expected, 

12 



given that in the present approximation the leading order term in Poo{x,z) corresponds to 
the asymptotic distribution of particles between two fiat plates. As a result, po does not 
contribute to the integral above. Now, using a Taylor expansion about A = for the domain 
of integration, and taking into account the periodicity of p„, i.e. L {dprjdx)dx = 0, we 
obtain the averge velocity to order A^, 



Uq = Pe; ui = 0; 



U2 = - 



dpi 



dx 



+ 



dpi 



z=l/2 



dx 



2 = -l/2 



gdx. 



(52) 



Let us note that, in agreement with the case of soft confinement, the first order vanishes, 
as expected given the symmetry of the problem. In fact, in both problems the geometry 
of the channel is not affected by a change A — ?■ —A and thus all the terms in the average 
velocity with odd powers of A are zero. Therefore, in order to calculate the leading order 
perturbation to the average velocity, U2, we first need to determine pi{x,z). 

Let us first calculate the perturbation to the marginal probability distribution for an 
arbitrary function g{x). Integrating the regular expansion for Pod{x, z) over the cross section 
and expanding the integral in a Taylor series about A, we obtain the following expression 
for the order A contribution to the marginal probability. 



1/2 

Pi{x) = 2g{x) + I pi{x, z) dz = 2g{x) + Si{x), 

-1/2 



(53) 



The governing equation for Si{x) is then obtained by integrating Eq. ( l50l) (for n = 1) with 
respect to z from —1/2 to 1/2 and using the corresponding no- flux boundary condition (see 
Appendix | 

0. (54) 



dx dx"^ dx 



In addition, Si{x) inherits the periodicity and normalization conditions from pi(x, z). The 
solution, after some manipulation, takes the form 

-/ e-'<g{OdC 

(55) 



si{x) 



-Pe 



-2g{x) + — 



2e 



Pex 



9{C)dC 

x+l 



Po 



dQ 



Therefore, the leading order perturbation to the marginal probability distribution is given 
by. 



pi(x) 



2e 



Pex 



-Pe 



- 1 



x+l 



dC 



(56) 
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This equation is analogous to Eq. (133|) . obtained in the case of soft confinement, where the 
perturbation to the width of the sohd channel, given by 2g{x), replaces the perturbation 
potential, /i(x). Moreover, taking into account that the unperturbed width of the channel 
considered here corresponds to Jq = 1, then the condition Ii{x) = 2g{x) is equivalent to 
the condition found in the case of narrow channels and discussed in Sec. Ill At that is 
/(x) = w{x). 

In order to obtain an analytical expression for the leading correction to the average ve- 
locity we consider the case g{x) = sin2'7rx. This choice oi g{x) satisfies the condition derived 
above, that is I{x) = w{x), for the special case of soft confinement with a perturbation 
given by Vi(x) = — 4sin27rx, discussed before in Sec. IIII All Clearly, the leading order 
perturbation to the marginal probability is then the same as in the case of soft confinement, 

2-7? sin 2'7rx — Pe cos 27rx 



Pi(x) =4tt 



(57) 



47r2 + Pe^ 

In Fig. [3] we compare the marginal probability density up to order A obtained from the 
asymptotic analysis with the results obtained from numerical simulations performed in both 
soft and solid cases for different Peclet numbers, where solid lines refer to the asymptotic 
results up to O(A^), and circles/stars refer to the results from the Brownian dynamics 
simulation in the soft/solid confinement. The agreement is excellent for Peclet numbers as 
large as Pe = 1. Let us note that the asymptotic analysis is valid for A ^ 1 and assumes 
that APe <^ 1 as well. For the cases shown in the figure this latter assumption is only valid 
if Pe ^ 10. The simulation results show that, as expected, the effect of soft and geometric 
confinement on the distribution of particles is not the same for large Peclet numbers, e. g. 
for large values of the driving force. 

The probabihty distribution pi{x,z) can be determined by proposing two base solutions 
of the form 

p^(a;,z) = /i(z)e2'^" and ^i(z)e-2™. (58) 

Then, substituting these general solutions into the governing equation and the equation 
for the no-flux boundary condition we obtain (see details in appendix ICj) . 

pi = 2R{Ai cosh(aiz)e^"^^), (59) 

where ]R(-) denotes the real part, and Ai and ai are the following complex numbers 

7r6 Pe 

Ai = — — -- and a? = e^(2mFe + An^). (60) 

aisinh(ai/2) ^ ^ ' ^ ^ 
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FIG. 3. Marginal probability density for different Peclet numbers. Solid lines refer to the asymp- 
totic results up to O(A^); circles/stars refer to the results from the Brownian dynamics simulation 
in the soft/solid confinement. Parameters used in the numerical simulation are A = 0.1 and e = 0.9. 



Finally, we calculate the average velocity at O(A^), 



U2 



— / 2 sin 27rx 







dpi 
dx 



dx = -in'^e Pe 



z=l/2 



coth(a/2) 



a 



(61) 



The corresponding effective mobility, normalized by the bulk mobility, is given by 



1^0 



Pe 



^ l-4AVe2: 



coth(a/2) 



a 



+ 0{X' 



(62) 



Comparing the equation above with Eq. |45] we see that the two mobilities are different. 
This is the first difference observed for the transport of Brownian particles between the 
cases of soft and geometric confinement. Moreover, even the limiting behavior at small 
Peclet numbers is, in general, different for the two types of channels. The limit of the 
effective mobility at zero Peclet number is. 



(v) 
lim ^^ = 1 - 2A27recoth(7re), 
Pc^o Pe ^ ^ 

and the asymptotic behavior at large Peclet numbers is 



(63) 



Pe 



1-A^ 



VvrPe 



(64) 
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FIG. 4. (Color online) Effective mobility of Brownian particles in solid and soft channels as a 
function of the Peclet number. The top plot corresponds to e=0.9 and the bottom one to e=0.5. 
In both cases we use A = 0.1. The solid (dashed) line corresponds to the effective mobility in 
solid (soft) channels accurate to order A'^. The open squares (circles) correspond to the results of 
Brownian dynamics simulations in the case of geometric (soft) confinement. 

Let us note, however, that the effective mobihties are equal in the hmit of narrow channels, 
as was already discussed in Sec. IIIAI In fact, taking the limit of e — ?■ we see that both 
effective mobilities tend to the same value, /ieff(e = 0) = 1 — 2A^. 

A comparison between the effective mobility in a solid channel and that in a channel 
created by a confining potential is presented in Fig. HJ The numerical results obtained from 
Brownian Dynamics simulations in both types of channels are also shown in the figure. The 
results from the numerical simulations agree very well with the leading order corrections 
given by Eqs. fH5|) and f l62|) . over a wide range of Peclet numbers. Moreover, in agreement 
with the asymptotic results obtained for large Peclet numbers, it is clear that the convergence 
of the effective mobility towards its bulk value is slower in the case of geometric confinement 
and, as a result, the mobility in solid channels is smaller at large Peclet numbers. On the 
other hand, the effective mobility in the case of geometric confinement is typically larger 
than that in the case of soft confinement for small values of the aspect ratio e. A detailed 
comparison of the mobilities in the cases of soft and geometric confinements as a function 
of both the aspect ratio e and the Peclet number is presented in Fig. O We can see that 
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FIG. 5. (Color online) Contour plot of the difference in the effective mobility between the cases of 
soft and geometric confinement, (/x|g — /x^)/;Uo. 

for aspect ratios above tc ~ 0.67, /Xgg > /x^ for all Peclet numbers. On the other hand, for 
smaller values of the aspect ratio, the difference in mobilities goes from negative (/i^g < /i^) 
at small Peclet numbers to positive as the Peclet number increases. For arbitrary but fixed 
values of the Peclet number the mobility difference goes from negative at small aspect ratios 
to positive as the aspect ratio increases. 

IV. CONCLUSION 



We have investigated the transport of Brownian particles confined by two-dimensional, 
weakly-corrugated channels. We considered the cases of soft confinement, in which the 
particles move in a periodic potential, and of geometric confinement between solid walls. In 
particular, we studied the effect that small variations in the width of the channels have on 
the average velocity of Brownian particles by means of regular perturbation methods. We 
showed that the leading order correction to the marginal probability distribution in the case 
of soft confinement is equivalent to that in the case of geometric confinement, provided that 
the (configuration) integral of the soft-confinement potential over the cross-section is equal to 
the channel width in the case of solid walls. The same condition was (previously) found in the 
case of narrow channels. We then considered the specific case of sinusoidal variations in the 
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channel width and calculated the leading order correction to the probability distribution and 
to the average velocity of the particles. We found that the reduction on the average velocity 
is different in the cases of soft and solid channels. Interestingly, whereas the reduction in the 
velocity is larger in the case of soft channels at small Peclet numbers and relatively narrow 
channels, the opposite is true at large Peclet numbers and for wide channels. The asymptotic 
behavior at large Peclet numbers is also different, with a much faster convergence to the 
asymptotic velocity in the case of soft confinement, with the perturbation decaying as 1/Pe^, 
compared to the case of geometric confinement in which the perturbation to the velocity 
vanishes as l/A/Pe. We also performed Brownian Dynamics simulations for the two types of 
channels and the numerical results agree well with the analytical results over a broad range 
of Peclet numbers. 
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Appendix A: Eigenvalue problem and eigenfunction expansion 

We first transform the linear differential equation for pi{x,z) given by Eq. ( I22l) into a 
self-adjoint linear differential equation for pi(x, z) = exp (— ^-|^) exp I -^ j pi{x, z), 



with 



52 g2 

C [pi{x, z)] = ( q{z) - e^— - — ) pi{x, z) = f{x, z), (Al) 



'W^'M|)%(if)^i^, 



fix,z) =pQ{z)ey 2 )e 



Poa;\ (Ml)." 



2 



2 



^-yw,^v/-^-v/ '''''' 



dx^ dz"^ V dz 



(A3) 



C is therefore a self-adjoint linear differential operator with an associated eigenvalue equation 
Cij) = Xijj that is separable into two linear ordinary differential equations. Specifically, 
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writting if) = XZ we obtain 

^X"{x) + I a - e^ (^^y I X{x) = 0, (A4) 

and the eigenvalue equation for C holds with A = a + /3. The first equation is a standard 
Sturm-Liouville problem in the interval (0, 1) and the eigenf unctions Xn{x) are sin(n7ra;) and 
cos(n7rx) with eigenvalues a„ = e^(n^'7r^ + Pe^/4). The spectrum of the second eigenvalue 
problem is also discrete, assuming that q{z) — > oo for z — )■ ±oo. We can therefore write the 
eigenfunctions and eigenvalues as Zm{z) and /?„, respectively. 

An interesting case is that of a quadratic potential, Vq{z) = ttz"^, confining the particles 
in the 2;-direction. In this case, the solutions to the second eigenvalue problem are Zm{z) = 
exp{—TTz'^/2) Hm{^/TTz), whcrc Hm{x) are the Hermite polynomials, and we can expand the 
general solution in known eigenfunctions, 

oo oo 

Pi{x,z) = ^ ^ [cn,msm{mix) + dn,mCos{mix)] e~^ Hm (V^z) , (A6) 

n=l m=l 



with 

dx / dz f{x,z) sm{mTx)e^2— H^(^y/^z), (A7) 



1 Z'+oo 



1 1 /"^ /-^^^ _^2 



"'" 2"n! j^^2^2 + ^) e2 + ^^ Jo 

dn,m = 7Z^,-/ J\ / d^ I dz f{x,z)cos{mTx)e^^ Hm{y/7Tz), (A8) 

^ '^- (n'^iT^ + ^je^ + 7imJo J-oo 



where f{x,z) is obtained from Eq. (IA3p . 

,d^Vi 



Pc X TT 2 



/(x, z) = e 2 e 2 



e^vrz^^^ + 27r (1 - 27rz^) Vi 



(A9) 



It is important to note that f{x, z) can be written as a combination of Zq{z) and Zi{z) and, 
therefore, the only non-zero coefficients are c„_o and c„_i. As a result, it is possible to write 
Pi{x,z) as a linear combination oi Pq{z) and po(-2)(27r2;^ — 1). 

Appendix B: Domain perturbation 

Here we apply the domain perturbation method^i^^ to transform the no-flux boundary 
condition at the curved wall defined hj z = h{x) = 1/2 + Xg{x) to an asymptotically 
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equivalent boundary condition applied aX z = 1/2. We start by substituting the regular 
perturbation expansion for the asymptotic probability density given by Eq. f l49|) into the 
no-flux boundary condition given by Eq. ([H]), 



Ae 



2^ 

dx 



Pe(po + ApO-^-A^ 
1 



+ ^ + A^ + A^%i + O(A^) = 0, 



dz 



dz 



dz 



(Bl) 



at z = — h Xg{x) 



Then, we approximate the functions p„(x, z) in the neighborhood of 2; = 1/2 with a Taylor 
series in powers of A, 



I _ \ ,^P^ 

Pn\z=l+\g(x) - Pn\z=\ + q^ 



1 (9Vn 



2 dz^ 



{\g{x)f + 0{\'). (B2) 



Substituting these regular expansions into the no-flux boundary condition and equating like 
powers of A we obtain the equivalent boundary conditions at z = 1/2 for each power of A, 

1 



at 0(1) 
at 0(A) 
at 0(A2) 



dpo ^ 
dz 

dx 
dx 



at z 



Pepo - 



dpo 
dx 



Pe pi + (? 



d'^pQ dpi ^ 1 



dpo 
dz 



9 



d'^Po dg dpo dpi 



, dp2 , d'^pi g^ d^po 



dxdz dx dz dx 



at z = -. 
2 



(B3) 



dz ' '' dz'^ 2 dz'^ 

Finally, we also need to transform the normalization condition in the entire channel to an 
asymptotically equivalent condition for a domain limited by the two parallel planes deflned 
by z = ±1/2. First, we substitute the regular perturbation expansion for the asymptotic 
probability distribution into the normalization condition. 



$(A) 



1 |.l/2+Ag(x) 

dx {po + \pi + X^p2 + 0{\^))dz = l. 

J-l/2-Xgix) 



(B4) 



Then, we expand the normalization integral $(A) in a Taylor series in the neighborhood of 
A = 0, 

(B5) 



$(A) = $(0) + A^ + -A^^ + O(A3). 



Calculating the expansion above explicitly, substituting it in the normalization condition, 
and equating like powers of A, we obtain the equivalent conditions on the distribution be- 
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tween the two parallel planes: 



at 0(1) : dx podz = 1; 

Jo J-l/2 

at 0(A): / ypo\^=i/2 + Po\;,=-i/2) 9dx+ dx pidz = 0; 



at 0(A^ 



^97 + '^^ 



2=1/2 



1/2 



z=-l/2 



(B6) 



(jfrfa; 



H H/2 

+ dx p2dz = 0. 

Jo i-l/2 



Appendix C: Leading order correction in the solid channel 

We seek a solution of Eq. (15 Op with the corresponding boundary and no-flux conditions to 
0(A). Given the linearity of the equation and the boundary and normalization conditions, 
we propose two base solutions of the form 



Pi{x,z) = fi{z)e''^''' and gi{z)e 



-2'Kix 



Substituting these proposed solutions into Eq. ( I50l) we obtain 

fi{z) = v4/icosh(ai2;), 
5-1(2;) = Agicosh{a^z), 

where ai is the complex number 



(CI) 

(C2) 
(C3) 

(C4) 



and ai is its complex conjugate. Since the probability distribution is a real value function 
we simplify the problem assuming Ai = Afi = Agi, 



Pi{x, z) = Ai cosh(a;i2;)e '"'^^ + Ai cosh(Q;iz)e 



-2nix 



(C5) 



This proposed solution clearly satisfles the periodicity condition in x, as well as the normal- 
ization condition. It is also symmetric about the 2;— axis, by construction. Therefore, the 
only remaining condition is the no-flux condition. Replacing then the proposed solution into 
Eq. (1B3P we determine the constant Ai, 

vre^Pe 



Ai 



ai sinh(ai/2) 



(C6) 
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